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Correlations between X-ray and radio spectral properties 
of accreting black holes 
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ABSTRACT 

We study correlations betwen X-ray spectral index, strength of Compton reflection, 
and X-ray and radio fluxes in accreting black holes (Seyferts and black-hole binaries). 
We critically evaluate the evidence for the correlation of the X-ray spectral index with 
the strength of Compton reflection and consider in detail statistical and systematic 
effects that can affect it. We study patterns of spectral variability (in particular, piv- 
oting of a power law spectrum) corresponding to the X-ray index-flux correlation. We 
also consider implications of the form of observed X-ray spectra and their variability 
for interpretation of the correlation between the radio and X-ray fluxes. Finally, we 
discuss accretion geometries that can account for the correlations and their overall 
theoretical interpretations. 

Key words: accretion, accretion discs - binaries: general - galaxies: Seyfert - radi- 
ation mechanisms: thermal - X-rays: galaxies - X-rays: stars. 



1 INTRODUCTION 

X-ray and soft 7-ray (hereafter X7) spectra from lumi- 
nous accreting black holes (hereafter BH), i.e., AGNs and 
BH binaries, commonly show a distinct component due to 
Compton reflection (Lightman & White 1988; Magdziarz 
& Zdziarski 1995) of the primary continuum from a cold 
medium (e.g.. Pounds et al. 1990; Nandra & Pounds 
1994; Magdziarz et al. 1998; Weaver, Krolik & Pier 1998; 
Zdziarski, Lubinski & Smith 1999, hereafter ZLS99; Done, 
Madejski & Zycki 2000; Eracleous, Sambruna & Mushotzky 
2000; Done et al. 1992; Gierlihski et al. 1997; Zdziarski et 
al. 1998; Zycki, Done, & Smith 1998; 1999; Gilfanov, Chura- 
zov & Revnivtsev 1999, 2000, hereafter GCROO; Revnivtsev, 
Gilfanov & Churazov 1999, 2001). A very interesting prop- 
erty of Compton reflection with a number of potential phys- 
ical implications is that its relative strength, ~ Q/2n, where 
Q is the solid angle of the cold reflector as seen from the 
hot plasma, correlates with some other spectral and timing 
properties of many sources (e.g., ZLS99; GCROO). 

Another correlation often found in both Seyferts and 
BH binaries is between X-ray spectral index and X-ray flux 
(e.g., Chiang et al. 2000; Done et al. 2000; Nowak, Wilms & 
Dove 2002; Zdziarski et al. 2002b, hereafter Z02; Lamer et al. 
2003; GUozzi, Sambruna & Eracleous 2003). The X-ray flux 
is also correlated with the level of radio emission in the hard 
states of BH binaries (Corbel et al. 2000; Gallo, Fender & 
Pooley 2003) . All these correlations appear to reflect funda- 



mental properties of BH accretion flows. We critically study 
the correlations, relationships between them, their theoreti- 
cal models, and the corresponding physical implications. In 
Appendix A, we present properties of spectral variability 
due to a power-law pivoting, which process is closely related 
to the flux-index correlation as well as it puts constraints on 
the interpretation of the radio-X-ray correlation. In spectral 
fits, we use xspec (Arnaud 1996). 



2 CORRELATION OF COMPTON 

REFLECTION WITH SPECTRAL INDEX 

2.1 Results of fits to data 

The first to find a correlation of Compton reflection with an- 
other spectral property were Ueda, Ebisawa & Done (1994), 
who found that Q correlates with the X-ray photon spectral 
index, P, in the BH candidate GX 339-4, albeit their result 
was based on only five observations (by Ginga). Later, 23 
Ginga observations of BH and neutron-star binaries were 
found to obey the same correlation (Zdziarski 1999). 

The reality of the correlation has been unambiguously 
confirmed in the RXTE data for the luminous BH binaries 
CygX-1, GX 339-4 and GS 1354-644 (Gilfanov et al. 1999; 
GCROO; Revnivtsev et al. 2001). Fig.Qpresents those RXTE 
results, as well as the Gmga results for 20 observations of 
CygX-1, GX 339-4 and Nova Muscae. These results were 
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Figure 1. Correlation between tiie strengtli of Compton reflec- 
tion and the X-ray spectral index in BH binaries. Small symbols 
with error bars and large ones without them correspond to ob- 
servations by RXTE (GCROO and references therein) and Ginga 
(Zdziarski 1999 and references therein), respectively. The solid 
curve corresponds to a model of ZLS99 with a central hot source 
surrounded by an overlapping cold disc. 



mostly obtained with spectra constrained to energies below 
a few tens of keV. However, in many instances, we also have 
at our disposal broad-band spectra extending to several hun- 
dred of keV, where almost the entire range of the reflection 
spectrum is found to fit the data very weU, together with 
thermal Comptonization (e.g., the hard states of Cyg X-1, 
Gierliriski et al. 1997, and GX 339-4, Zdziarski et al. 1998), 
as illustrated in Fig.|5| 

The Q,-V correlation is also seen in Fourier-resolved 
spectra, i.e., corresponding to variability in a given range of 
Fourier frequencies (Revnivtsev et al. 1999, 2001; Gilfanov 
et al. 1999). Very importantly, the strength of reflection was 
also found to correlate with the low-frequency QPO centroid 
frequency and with the degree of the relativistic smearing of 
the Fe Ka line (associated with reflection, Zycki & Czerny 
1994), see flg. 4 in Gilfanov et al. (1999) figs. 1.4 and 1.6 in 
GCROO, fig. 2 in Revnivtsev et al. (2001), and Gilfanov et 
al., in preparation. 

We note that the above results concern BH binaries in 
their luminous (and mostly hard) states. Existing data are 
insufficient to constrain reflection in quiescent states of BH 
binaries. However, given the prevalent theoretical interpre- 
tation of the correlation as due to feedback between a cold 
accretion disc and hot plasma (Section 15. 11 . we do not ex- 
pect it to be present in quiescence as the disc is then cut 
off at a large radius (e.g., Narayan & Yi 1995) and Comp- 
tonization of synchrotron photons usually dominates X-rays 
(e.g., Wardziriski & Zdziarski 2000). Analogously, we do not 
expect an Q.-T correlation in low-luminosity AGNs. 

The picture of the correlation disappearing with the de- 
creasing luminosity is indeed conflrmed in the BH transient 
GS 2023-1-338. The results of Zycki et al. (1999) show an Q,-V 
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Figure 2. A broad-band ~ 1-10"' keV spectrum of Cyg X-1 
(Gierliriski et al. 1997) fitted by thermal Comptonization and 
Compton reflection (black solid curve). The two components sep- 
arately are shown by the short dashes and dots, respectively. At 
low energies, the spectrum is absorbed in interstellar and circum- 
stellar media. The long dashes and dot/dashes show (unabsorbed) 
nonthermal synchrotron spectra emitted by power-law electrons 
with p = 2.5 and with an exponential and sharp high-energy cut- 
off, respectively (see Section l4l. 



correlation in the initial luminous hard state, but then the 
spectrum substantially softened at an approximately con- 
stant f2 when the luminosity decreased by a factor > 10. 
This effect is most likely due to the onset of the dominance 
of Comptonization of synchrotron photons with decreasing 
luminosity (Wardzinski & Zdziarski 2000). Similarly, reflec- 
tion is weak in the transient XTE J1118-f480 (Frontera et 
al. 2001b; Miller et al. 2002), which maximum luminosity 
was only about ~ 10^^ of the Eddington one^. 

The first to report an Q.-V correlation in an AGN were 
Magdziarz et al. (1998), in the Seyfert-1 galaxy NGC 5548. 
Then, ZLS99 showed the presence of a strong Q.-T correla- 
tion in 47 Ginga observations of 23 AGNs (mostly Seyfert 
Is and a few intermediate-type AGNs). Matt (2001) pre- 
sented a compilation of BeppoSAX results for Seyfert Is, 
and found that Q is correlated with F in the full sample at 
the confidence level of > 0.999, and ~ 0.99 if one object 
in which the measurement of the continuum was probably 
affected by a soft excess is not taken into account. A similar 
sample was later studied by Perola et al. (2002, hereafter 
P02b). Zdziarski & Grandi (2001, hereafter ZGOl) showed 
that Compton reflection in broad-line radio galaxies (also 
called radio-loud Seyfert Is), albeit weaker on average than 
in radio-quiet AGNs, is still consistent with the same cor- 
relation. Nandra et al. (2000, hereafter NOO) found an Q-T 
correlation in multiple RXTE observations of the Seyfert 
galaxy NGC 7469. Papadakis et al. (2002, hereafter P02a) 



^ An additional complicating factor in measuring reflection in 
this halo system is its likely overall low metallicity, which possi- 
bility was considered by Frontera et al. (2001b), but not by Miller 
et al. (2002), who only considered the case of a low Fe abundance 
at all other abundances kept solar. Thus, the actual value of Q/2n 
in this object may still be ~ 0.2 found by the former authors. 
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Figure 3. Correlation between the strength of Compton reflec- 
tion and the X-ray spectral index (in the ~ 3-20 keV range, see 
Section 12.41 and eq. 0) in AGNs. Open and filled symbols corre- 
spond to radio-quiet AGNs and active radio galaxies, respectively. 
Blue circles, green triangles, and red squares correspond to ob- 
servations by Ginga, RXTE, and BeppoSAX, respectively. Error 
bars/contours are not shown for clarity, see Fig.|3J. for Ginga er- 
ror contours. References for Ginga: Lubiriski & Zdziarski (2001), 
Wozniak et al. (1998); RXTE: Weaver et al. (1998), Lee et al. 
(1998), Done et al. (2000), Chiang et al. (2000), P02a, Eracloous 
et al. (2000), this work (NGC 7469 in Section IT^: BeppoSAX: 
P02b, Orr et al. (2001), ZGOl, Grandi et al. (2001). 



have found this correlation in the average properties of four 
Seyferts. Fig. |21 summarizes those (and some other) results. 

It is, however, of importance to consider possible effects 
that may lead to spurious fl-T correlations. We consider be- 
low in detail statistical and systematic effects. 



2.2 Statistical effects 

The main statistical effect is due to the fitted strength of 
reflection and spectral index for a single observation being 
correlated to certain degree, resulting in a skewness of their 
joint error contour. Any data have limited statistics, and 
then the same intrinsic ~ 2-20 keV spectrum can be fit- 
ted within some confidence limit with either a somewhat 
harder index and less refiection or a softer index and more 
refiection. This effect is more important for AGNs than for 
BH binaries, as the former have usually significantly lower 
statistics. 

As usual with this type of effects, its importance de- 
creases with increasing statistics. As shown in Fig. the 
extent of typical error contours for Seyferts is much less than 
the global extent of the correlation, and the inclination of 
an individual contour is also significantly steeper from that 
of the correlation. Thus, it is already very unlikely that sta- 
tistical effects alone can account for this distribution. 

Then, Vaughan & Edelson (2001, hereafter VEOl) con- 
sidered 120 simulated RXTE spectra based on observations 
of MCG -6-30-15 with short exposure times of ~ 1.4-3.7 
ks and the top layer of the RXTE PCUs 0-2. They used a 
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Figure 4. (a) Error contours {la for 2 parameters) for the cor- 
relation between Q and F in the Ginga data of ZLS99 for AGNs. 
For clarity of display, only contours with the vertical extent of 
A(S7/27r) < 1 are shown, and alternating colors and line widths 
are used. The panels below show simulations, see text for details, 
(b) Best fits and 2 representative contours for simulated data with 
a uniform distribution of F and fixed Q. (c) Fixed Q for 6 values of 
F. (d) Random Q for a uniform distribution of F. Clearly, all sim- 
ulations results are very different from the observed distribution. 



power- law model with F in a relatively narrow range, ~ 1.9- 
2.15 for 90 per cent of the simulated spectra, and constant 
reflection at ^1/2tt = 1.42. They fltted the simulated spec- 
tra with r and fl/2n as free parameters. The fltted values 
of (F, Q/2n) formed a skewed elongated contour extending 
from ~ (1.75, 0.5) to ~ (2.3, 3). 
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Given the assumed narrow range of F, the (F, Q./2ti) 
contour obtained by them is equivalent (except for some 
spread in the values of F and the exposure time) to the sta- 
tistical error contour for two parameters. Indeed, the typical 
error bars provided by VEOl are comparable to the extent of 
their contour. However, based on the extent of this contour, 
VEOl stated that the above statistical effect casted serious 
doubts on the claim of the correlation by ZLS99. 

However, most of the error contours of ZLS99 (Fig. 
and of those corresponding to RXTE and BeppoSAX re- 
suhs of Fig. 01 are much smaller than the contour obtained 
by VEOl due to much better statistisc. In particular, the 
typical statistics in the Gtnga spectra is an order of mag- 
nitude better than those of the 1-orbit simulated spectra 
of VEOl. Also, the shape of the error contour of VEOl is 
significantly difi'erent from those observed (Fig.j^J. Already 
based on the above, the suggestion of VEOl is unlikely to be 
correct. 

In order to quantitatively study this issue, we have gen- 
erated 208 simulated RXTE spectra (assuming a 0.5 per 
cent systematic error, which we estimated from our fits to 
Crab data) with the range of F of 1.90-2.15 (similar to that 
of VEOl) and the constant f2/27r = 1. In the simulations, we 
have assumed a source with the constant 2-10 keV fiux of 
3.5 (flux unit is 10~^^ erg cm~^ s~^ hereafter in this para- 
graph) and the exposure of 50 ks for PCUs 0-2 and 30 ks 
for PCUs 3-4, similar to that of observations of NGC 7469, 
see Section r2.5l below. The product of the flux and exposure 
roughly corresponds to typical RXTE measurements shown 
in Fig. 13 For example, NGC 5548 has the flux in the range 
of ~ 5-10, exposures of ~ 20-40 ks (Chiang et al. 2000; see 
also Section below). IC 4329A has the flux of ~ 15, ex- 
posure of ~ 70 ks (Done et al. 2000), and 3C 120 has the 
flux of ~ 6, exposure of ~ 60 ks (Eracleous et al. 2000). 
The results are shown in Fig. We clearly see that the 
obtained distribution is very difi'erent from those in Fig. |H] 
and|^. It is also different from the contour of VEOl due to 
the better statistics. 

In Fig. ^h, we also show two representative error con- 
tours. We see that although they correspond to only la, 
they do provide a good measure of the spread of almost 
all the simulated points. This overestimate of statistical er- 
rrors by xspec is due to the corresponding overestimate of 
background errors (added in quadrature to the total ones) 
discussed by NOO and VEOl. In order to account for this ef- 
fect, NOO did not include background errors in their data at 
all, and VEOl performed simulations to estimate the actual 
errors. In the case of the Ginga data of ZLS99, the stan- 
dard XSPEC procedure was followed. Thus, the error con- 
tours on Fig.|l^ are somewhat overestimated. However, this 
is a conservative approach, which can only reduce the mea- 
sured strength of an actual correlation but will not lead to 
an appearance of a spurious one. Therefore, we flnd the sug- 
gestion of Edelson & Vaughan (2000) (based on the results 
of VEOl) that the correlation of ZLS99 could be spurious 
due to inadequate modeling of errors to be incorrect. 

Results of simulations with a larger range of the input 
F (and other assumptions similar to those above except for 
the systematic error being now 0.01) are shown in Fig. |1J;. 
These results also show how the fitted values of F and Q./2-n 
correlate with the ones assumed. We see no systematic ef- 
fects here apart from the dispersion of the fitted value of 



Q./2-n increasing with the increasing F (due to decreasing 
photon statistics at high energies). 

To further study possible contributions from statistical 
effects, we have also performed analogous simulations (us- 
ing the exposure time, 5.4 ks, and response as for a Ginga 
observation of NGC 3516^) assuming Q,/2-r to be random in 
the 0-1.5 range. Results for 100 simulated spectra are shown 
in Fig.|lji. We see that the obtained distribution is still very 
different from that of Fig. UK- 

Our results confirm the finding of ZLS99 that the prob- 
ability that a spectrum with (F, ^/2n) measured to be at 
the high end of the extent of the correlation shown in Fig. 
2Ji would correspond in reality to a point at a low end of the 
correlation is extremely low. ZLS99 noted that (F, 57/27r) 
contours are skewed and elongated (see Fig. 2t)i s^nd thus 
used a statistical method taking into account the correlated 
errors while looking for the actual functional dependences 
between parameters. Based on comparison of the resulting 
values of (Bevington & Robinson 1992), ZLS99 deter- 
mined the probability of the correlation appearing by chance 
as ~ lO"'^". The exact value may be somewhat different be- 
cause the model with an allowed ^(F) dependence yields 
/v ~ 2 (which reflects the intrinsic spread of the data). 
Still, the probability of reducing x^ from 312 (assuming con- 
stant f2) to 115 (allowing for a phenomenological S7[F] de- 
pendence as a power law) by chance is in any case ^ 1. 

2.3 Intrinsic spread and small data samples 

We see in Figs. 13 and 2t that apart from the correlation 
between f2 and F, there is also a clear intrinsic spread of 
the values of those parameters in the observed F-f2 space. 
This spread is caused by varying model assumptions used in 
the references on which Figs. Q |H] are based, differences in 
calibration of different instruments, and, last but not least, 
fluctuations due to some physical effects (e.g., orientation) 
superimposed on the overall n(F) dependence. The last ef- 
fect is rather common in astrophysics, in which it is very 
rare that any two quantities are correlated without any in- 
fluence of other parameters of the system. In the case of BH 
binaries, we see a full width of A(fi/27r) ~ 0.2 or so for a 
given F, whereas it is ~ 0.5 for Seyferts. 

Thus, a finding of no fi-F correlation in a small sample 
cannot be taken as a proof of the lack of it, even if the sta- 
tistical quality of data is very good. Furthermore, a sample 
with poor statistics cannot be considered a basis of a proof 
for either presence or absence of the correlation. For ex- 
ample, the data points for GX 339-4 by Wilms et al. (1999) 
have so large error bars and so limited range of F and Q. that 

^ The Ginga spectra available to us are already background- 
subtracted, and their errors have been calculated including the 
uncertainty of the background. On the other hand, simulations 
give errors neglecting the effect of the background. For this par- 
ticular source, this results in a decrease of the errors roughly 
corresponding to an increase of the exposure by a factor of ~ 2. 
Thus, this set of simulations corresponds to the actual exposure 
of ^ 10 ks, which is still a half of the actual average exposure in 
the data of ZLS99 of 20 ks. On the other hand, the effect of back- 
ground errors in XSPEC is generally overestimated (NOO, VEOl), 
and thus the present simulations may better represent the actual 
uncertainties than those in Figs. Hb-c. 



© 2003 RAS, MNRAS OOO.nTOl 



Correlations between X-ray and radio spectral properties of accreting black holes 5 



no statement about the correlation can be made at all. In 
the case of that source, the actual presence of a strong corre- 
lation (see Fig. has later been established by Revnivtsev 
et al. (2001) and Nowak et al. (2002). Then, P02b claimed 
the fi-F correlation to be relatively weak in their sample 
of Seyferts. However, the error bars in their data are large 
enough for their data to permit the actual presence of a 
strong correlation. The same holds for the similar data set 
of Matt (2001), who raised the issue of the influence of sta- 
tistical effects on the ^(F) correlation found by him in the 
BeppoSAX data but presented no calculations to test it. 

2.4 Systematic effects 

Various systematic effects might affect the best fit values of 
the Q and F and can lead to biased estimates and to con- 
tradiction between results obtained by different instruments 
and using different spectral models. In certain circumstances 
these effects could potentially lead to appearance of spurious 
corelations between parameters. These effects are considered 
below. 

NOO noted that a background subtraction with a sys- 
tematic relative error with a specific power-law dependence 
on energy (reaching 0.03 at 20 keV) may produce a spu- 
rious Q-F correlation in their RXTE data for NGC 7469. 
However, it is easy to see that this would be the case only 
if the dominant variability pattern corresponded to a pivot 
point at low energies, < 10 keV. Apart from the specific 
case of NGC 7469, we show in Section |3 that the typical 
pivot energy in Seyferts is S> 10 keV, which, in turn, would 
produce an Q,-T anticorrelation. Thus, this explanation can- 
not be general. Also, in Section [2. 51 below, we reanalyze the 
data of NOO with an updated PGA background model and 
recover the same correlation. Furthermore, we consider it ex- 
tremely unlikely that background would be undersubtracted 
in the same way in AGNs observed by Ginga, BeppoSAX, 
and RXTE. 

Then, Weaver et al. (1998), and later P02b and Malzac 
& Petrucci (2002, hereafter MP02), noted that a correction 
needs to be made to account for the high-energy cutoff in 
the spectra. Namely, the lower the cutoff the fewer incident 
photons are available for reflection, and then the fitted value 
of n/27r increases. A related important effect not noticed by 
P02b is the downward curvature of the incident model spec- 
trum (present in the e-folded power law but not necessarily 
in Comptonization models) being compensated by an in- 
crease in the reflection strength. These two effects help to 
explain the offset between the Ginga and BeppoSAX results 
at low values of reflection noticed by Lubihski & Zdziarski 
(2001), and seen in Fig. 13 On the other hand, P02b notice 
that the resulting correction to cannot by itself remove 
the presence of the correlation. 

We also point out here the existence of an important 
correction to the X-ray index not taken into account by 
P02b. It is the difference between the asymptotic low-energy 
index, Ff, of an e-folded power law. 



AF = 



E2 - El 
lnE2/Ei' 



E7' ~ ^ 



9_ 



(2) 



— ccE ^exp 



E_ 
E'i 



(1) 



(where A'^ is the photon number) , and the actual X-ray index 
between two energies, Ei and E2. The latter equals F = 
Ff + AF, where 



and the second equality corresponds to i5i = 3 keV (limiting 
the range with possible dominance of a soft X-ray excess or 
effects of ionized absorption) and E2 — 20 keV, used by 
ZLS99. At E{ = 400 keV assumed by them, AF ~ 0.02, 
which is negligibly small. Thus, the correlation presented in 
ZLS99 concerns the actual hard X-ray index. 

On the other hand, this model fitted to BeppoSAX data 
on Seyferts by P02b (see also Matt 2001) yield generally 
lower values of Ef. Then, the fitted values of i}/2n increase 
somewhat with respect to those corresponding to fits with 
no cutoff. For example, for E{ — 160 keV, the fitted value of 
Q increases by ~ 50 per cent (P02b). However, as pointed 
out above, the value of the index needs to be corrected as 
well, with AF = 0.06. In fact, the two points P02b found 
to disagree with the correlation of ZLS99 both require sub- 
stantial corrections in F. In particular, E{ = 67 keV fitted by 
P02b for Mrk 509 yields a rather large AF ~ 0.13, increas- 
ing the 3-20 keV index from Ff = 1.58 (P02b) to F ~ 1.71. 
Taking into account both corrections, to Q and F, makes 
the BeppoSAX results (Matt 2001; P02b) quite compatible 
with those from Gmga of ZLS99 (see Fig. |3J . 

Furthermore, we point out that an e-folded power law 
is actually a very bad model for thermal Comptonization, 
which physical model fits well the spectra of Seyferts at high 
energies (e.g., Zdziarski, Poutanen & Johnson 2000). For 
spectral parameters characteristic to both Seyferts and BH 
binaries, thermal Comptonization has a much sharper cutoff 
(following an extended, power- law like, part of the spectrum) 
than that of the e-folded power law. This is illustrated in 
Figs. 1^, b for two cases approximately spanning the range 
of parameters of P02b. This further demonstrates that Ff of 
the e-folded power law model has no direct physical mean- 
ing. The inadequacy of the e-folded power law model was 
also noted by MP02, who have found that the sharper cut- 
offs of the thermal Comptonization spectra lead to the in- 
crease of the value of fl being substantially less than that 
for the e-folded power law. 

Concluding, results of fits with e-folded power laws need 
to be treated with caution. We also note that Zdziarski et al. 
(1995) and Gondek et al. (1996) obtained values of Ef > 500 
keV for the average Seyfert spectra using data from the 
CG-RO/OSSE detector (which affords coverage to signifi- 
cantly higher energies than BeppoSAX), i.e., substantially 
higher than those from BeppoSAX. 

In general, there are many other potential systematic 
effects that can affect the absolute determination of r2/27r, 
e.g., due to an approximate treatment of ionization of the 
reflecting medium. Then, in order to test the reality of the 
correlation, we should ask whether a given set of fits prop- 
erly ranks the spectra by the reflection strength. This issue 
was addressed by Gilfanov et al. (1999, see their fig. 7) and 
Revnivtsev et al. (2001, see their fig. 3) in the case of Gyg 
X-1 and GX 339-4, respectively. They have shown that the 
ratio of the count spectrum fitted with a higher value of 
Q/2n to that fitted with a lower value itself has the shape 
characteristic to Gompton reflection. This confirms that fits 
using relatively simple reflection models do provide a correct 
ranking in the values of 

We would also like to comment on the use of a nor- 
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Figure 5. Spectra from isotropic thermal Comptonization com- 
pared to e-folded power laws. Blackbody photons with fcT^i-, = 10 
eV distributed uniformly in a sphere are Comptonized by ther- 
mal electrons with the temperature, fcT, and the radial Thomson 
depth, T. The dotted, solid and dashed curves correspond to the 
distribution of the input blackbody, Comptonized spectrum cal- 
culated by the method of Poutanen & Svensson (1996), and a 
characteristic e-folded power law with the low-energy index of 
Ff and the folding energy of Ef. The filled squares with vertical 
error bars correspond to Monte Carlo results using the code of 
Gierhnski (2000). The parameters are: (a) kT = 100 keV, r = 1, 
Tf = 1.85, E{ = 300 keV; (b) kT = 25 keV, t = 4, Tf = 1.65, 
E{ = 150 keV. We see that in neither case any e-folded power law 
can approximate well the actual spectrum. 



malization of the reflection spectrum of Nandra & Pounds 
(1994) and NOO. They characterized the strength of reflec- 
tion by the 1-keV flux, Arof, of a power law with the same F 
as the incident one but normalized to n/2n = 1. Although 
equivalent to using the measured Q/2tt, this definition may 
lead to a spurious systematic correlation between Aj-d and F. 
To illustrate it, let us consider fl/2'K = constant. A common 
variability pattern of the X-ray power law spectra of both 
AGNs and BH binaries is variable F with an approximately 
constant pivot energy, Ep, 



dE [ E^ 



(3) 



where C is the normalization at Ep, and usually £p 2> 1 
keV (see Section |21 below). Then, 



oc {Ep/lkeVf, 



(4) 



i.e., Aref is strongly positively correlated with F even if the 
actual reflection solid angle remains constant. Thus, we ad- 
vise against using this measure of reflection. (We note, how- 
ever, that the specific data for NGC 7469 of NOO show F 
uncorrelated with the X-ray fiux, in which case the above 



effect is not critical and probably does not affect substan- 
tially results of that paper.) 

2.5 Single objects vs. samples 

Another issue to bear in mind when considering evidence for 
and against the Q-F correlation are the distinctions between 
different classes of sources and between samples of sources 
and repeated observations of a single source. The evidence 
for the presence of this correlation in Seyferts was given by 
ZLS99 mostly for broad-line Seyfert Is as a class. Repeated 
observations of single Seyferts give mixed results, sometimes 
showing the correlation in a given source (e.g., Magdziarz et 
al. 1998 for NGC 5548; NOO for NGC 7469, also see below), 
and sometimes not, e.g., in IC 4329A (Done et al. 2000) and 
some other sources (P02b). 

Fig. shows Q-T contours for RXTE observations 
the Seyfert 1 NGC 7469 (the inclination was assumed to 
be i = 30°). The data are the same as those of NOO, but 
they have been reextracted by us using the LHEASOFT 5.2 
version of the PCA response matrix and the background 
model. We used PCUs 0-4 whenever available and grouped 
the data in 4-day segments, and included a 0.5 per cent sys- 
tematic error. This yields the exposure times of ~ 50-70 and 
~ 30-50 ks for the PCUs 0-2 and 3-4, respectively. The con- 
tours clearly show a correlation (which would become even 
stronger after the reduction of the statistical errors discussed 
by NOO and VEOl, see also Section [2.21 . In order to test its 
reality, we have generated > 100 simulated spectra (using 
the same exposures and fluxes as for the actual data), a se- 
lection of them shown in Fig. |^ for two assumed values of 
F. In all cases, we flnd the statistical n(F) dependence to 
be much steeper than the observed one. Thus, we confirm 
the corresponding conclusion of NOO, who also ruled out the 
origin of the correlation in NGC 7469 from statistical effects 
(see their appendix). 

Fig. shows the results for NGC 5548 from Ginga 
by ZLS99 (contours) and from RXTE by Chiang & Blaes 
(2003) (error bars). Note good agreement between both sets 
of measurements after the update of the RXTE instrumen- 
tal response by Chiang & Blaes (2003) with respect to the 
original result of Chiang et al. (2000), who claimed a dis- 
agreement with ZLS99. 

In the case of BH binaries, the correlation is seen both 
in a number of individual objects (Cyg X-1, GX 339-4, Nova 
Muscae in the hard state) as well as in those sources consid- 
ered together (Fig.0. On the other hand, it is certain that 
there are some BH binaries that do not obey the 57-F corre- 
lation. In particular, a disapearance of the correlation with 
a decreasing Eddington ratio was pointed out in Section r2.1l 

Although the data for Cyg X-1 in the soft state lie on 
the extrapolation of the dependence for the hard state for 
this and other sources (Fig. see also fig. 5 of Gilfanov 
et al. 1999), the existing data appear insufficient to conclu- 
sively show the presence or absence of the Q.-T correlation 
within the soft state. On the other hand. Ran & Greiner 
(2003) claimed the presence of an Q.-T correlation in rela- 
tively soft states of GRS 1915-1-105. Then, Ballantyne, Iwa- 
sawa & Fabian (2001) did not find such a correlation in their 
sample of 5 narrow-line Seyfert Is, which class of object is 
likely the extragalactic counterpart of the soft state of BH 
binaries (Pounds, Done & Osborne 1995). Similarly, Lamer 
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Figure 6. f! vs. F in two Seyferts. (a) The 1-a contours for NGC 
7469 obtained by us for the RXTE data (NOO). (b) Simulations 
for NGC 7469. We confirm the conclusion of NOO that the extent 
of the correlation cannot be explained by statistical effects, (c) 
Contours for NGC 5548 from Ginga (ZLS99), and error bars from 
RXTE (Chiang & Blaes 2003). 



et al. (2003) found no Q-T correlation in the narrow-line 
Seyfert NGC 4051. The lack of the correlation in this object 
may be due to the domination of the emission process by 
non-thermal electrons (see Section fS. It . which possibility in 
soft-state sources was pointed out by ZLS99. 



3 X-RAY INDEX-FLUX CORRELATIONS 

Another common correlation concerning spectra of accret- 
ing BHs is that between the X-ray spectral index, F, and 
the X-ray flux, F. A common situation in the hard state 
of BH binaries and Seyfert Is is F showing an increasing 
trend (although with often significant non-statistical scat- 
ter) with F in an X-ray energy range, E\-E2- Examples of 
this behaviour is shown in Fig.|7|for the Seyferts NGC 5548 
(Chiang & Blaes 2003) and NGC 4051 (Lamer et al. 2003). 
The T-F dependence can often be fitted by a power law (e.g., 
Chiang et al. 2000; Done et al. 2000), also for the data shown 
in Fig.|7] The linear dependences of r(logF) in Fig.Qyield 
/v < 1, i.e., the departures from the linear correlation are 
compatible with being statistical only. On the other hand, 
the assumption of a constant F results in /v = 25/6 and 
138/6, respectively, and the F-test (Bevington & Robinson 
1992) yields the probability that the fit improvement is by 
chance of ~ lO"'^ and ~ 10~^, respectively. Thus, the V-F 
correlations are highly significant statistically, which conclu- 
sion is also confirmed by application of the Spearman and 
Kendall rank correlation tests. Similar power-law correla- 
tions are found, e.g., in the hard state of Cyg X-1 on long 
time scales (Z02). 

If the index increases from Fh to Fs while the energy 
fiux increases from to Fs, the two power laws intersect 
at the pivot energy. 



Fs /2-F, 



Fh V2-Fhy (F2/F1) 



(F2/F1)' 



i/(rs-rh) 



with the following substitution if either F equals 2, 



{E2/E^f 



1 



In 



E2 

'eI 



(5) 



(6) 



If Fs/Fh represents the ratio of the photon, rather than en- 
ergy, fluxes, then 2 — F above should be replaced everywhere 
by 1 — F. The resulting variable power law is given by equa- 
tion JHJ. A treatment of the effect of pivoting on average 
spectra and spectral variability is given in Appendix A. 

In the cases of the data for NGC 5548 and NGC 4051, 
Ep ~ 180 keV and 140 keV, respectively. As shown in Ap- 
pendix A, the assumption of a constant Ep results in an 
almost linear dependence between F and the logarithm of 
the flux over some energy band. Indeed, the E1-E2 energy 
flux, given by 



Fe,^ 



CEl^ 



El 



2 - F 



(7) 



and plotted in dashed lines in Fig. |7| almost coincide with 
the best-flt linear dependences. 

Table 1 gives a number of examples of the pivot en- 
ergy in Seyferts, as well as in Cyg X-1 and Cyg X-3. Table 
1 is based mostly on figures in published papers, and thus 
the numbers there are approximate only. Still, they unam- 
biguously show that many Seyferts have the pivot energy 
at ^ 10 keV. A similar pattern is shown by MCG -5-23-16 
(Zdziarski, Johnson & Magdziarz 1996). Although variable 
X-ray absorption may affect the correlation for NGC 4151 
found by Yaqoob et al. (1993), the weak variability found 
above 50 keV (Zdziarski et al. 1996; Johnson et al. 1997) is 
consistent with the pivot at high energies. 

Table 1 gives two examples of AGNs with the pivot at 
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Table 1. Pivoting in accreting black holes. 



Object 






r. 






Reference 




[kcV] 








[kcV] 




Cyg X-1 (hard state) 


3-12 


1.4 


2.0 


3 


40° 


Z02 


Cyg X-3 (hard state) 


3-100 


NA 


NA 


NA 


~ 20^= 


McCoUough et al. {1999b) 


3C 120 


0.3-2 


1.70 


2.03 


2d 


5 


ZGOl 


3C 390.3 


1-10 


1.70 


1.90 


1.8^ 


90 


Wozniak ct al. (1998) 


3C 390.3 


2-10 


1.61 


1.79 


2.1 


300 


Gliozzi ct al. (2003) 


IC 4329A 


2-10 


1.90 


2.07 


1.9 


200 


Done et al. (2000) 


IC 4329A 


2-10 


1.75 


1.95 


2.2 


240 


Madejski et al. (2001) 




2-10 


1.8 


2.2 


3 3d 


80 


VEOl 


MCG -6-30-15 


3-10 


1.9 


2.2 


3.8/ 


410 


P02a 


Mrk 766 


1-10 


1.64 


2.01 


1.5= 


10 


Leighly et al. (1996) 


NGC 3227 


2-10 


1.49 


1.75 


1.4 


20 


Ptak et al. (1994) 


NGC 3516 


2-10 


1.63 


1.69 


1.37 


910 


Chiang (2002) 


NGC 4051 


3-10 


1.3 


2.6 


26^' 


60 


P02a 


NGC 4051 


2-10 


1.60 


2.33 


12.3 


140 


Lamer et al. (2003) 


NGC 4151 


2-10 


1.4 


1.7 


6 


2000 


Yaqoob et al. (1993) 


NGC 5506 


3-10 


1.9 


2.1 


4/ 


5000 


P02a 


NGC 5548 


2-10 


1.65 


1.80 


1.73 


180 


Chiang & Blaes (2003) 


NGC 5548 


3-10 


1.8 


2.0 


3.2-f 


1700 


P02a 



Energy flux ratio except when noted. 
^ Also observed in the broad- band, 1.5-300 keV, variability (Z02). 
Implied by the 3-100 keV variability. 

Count rate ratio assumed here to represent photon flux ratio. 
^ Calculated using the 1-keV normalization. 

/ The (3— 5)+(7— 10) keV count rate ratio assumed here to approximate the 3—10 keV photon flux 
ratio. 



< 10 keV, 3C 120 and Mrk 766. In the cases of 3C 120 as 
well as Cyg X-1 in the hard state (where the seed photon 
energies are much higher than those in AGNs, which results 
in the pivot energy higher than that in 3C 120, see Z02), 
ZGOl and Z02 have shown that the spectral variability is 
consistent with the bolometric luminosity of the high-energy 
source being constant. This type of behaviour may also take 
place in the Seyfert IH 0419-577, where Ep ~ 7 keV (see 
fig. 3 in Page et al. 2002). However, Table 1 shows that this 
type of behaviour is certainly not the rule. The pivot at high 
energies implies the bolometric luminosity increasing with 
the softening of the source, as illustrated in Fig. |H| This 
is also confirmed by the 1-300 keV luminosity estimated 
by P02a for MCG -6-30-15, NGC 5506 and NGC 5548 to 
increase by a factor of ~ 3 for V increasing within the ranges 
given in Table 1. A similar increase is shown for NGC 5548 
in fig. 6a of Magdziarz et al. (1998). The results of Nowak et 
al. (2002, see their table 1 and fig. 3) imply the bolometric 
luminosity is increasing with the increasing 3-9 keV fiux in 
the hard state of the BH binary GX 339-4. Some increase of 
the bolometric luminosity with increasing F is also found in 
Cyg X-1 as the source gets close to a transition to the soft 
state (Z02). 

In the case of MCG -6-30-15, Shih, Iwasawa & Fabian 
(2002) have found a V-F correlation in ASCA data to appear 
to saturate above certain count rate. They have interpreted 
this behaviour as due to a superposition of two power law 
components with constant indices, and the normalization of 
the softer power law was allowed to vary. This model also 
allows to explain an approximate constancy of the Fe Ka 



line fiux with varying continuum fiux in those data. Thus, 
this is an attractive model for MCG -6-30-15. On the other 
hand, data for many objects clearly show correlated hard- 
ness changes over rather broad bands (see the references 
above, e.g., P02a), which also rules out the above interpre- 
tation as general. Furthermore, the apparent saturation at 
a high F seen by Shih et al. (2002) (and Merloni & Fabian 
2001, hereafter MFOl) might be an artifact of plotting F 
against the linear fiux^ rather than log F. As discussed above 
(see also Appendix A), variability with a constant pivot en- 
ergy results in a linear dependence between F and log_F, 
not F itself. Indeed, the nonlinearity and saturation of F in 
NGC 4051 claimed by Lamer et al. (2003) is specific to us- 



^ Often, the ratio of instrument counts in two bands (hardness 
ration) is plotted as a function of a count rate. The advantages 
of this choice are independence of possible revisions of the instru- 
ment response (giving the conversion between counts and pho- 
tons), which revisions are relatively common in X-ray astronomy, 
and no need to assume a spectral model. Another choice is to 
plot the fitted spectral index against the fitted energy flux. Al- 
though both quantities depend on the response, the advantage of 
this choice is an ease of comparison with physical models. On the 
other hand, VEOl, MFOl and Shih et al. (2002) show hybrid plots 
with the fitted F against the instrumental count rate. The last two 
papers show a comparison of the results with a physical model by 
MFOl assuming the X-ray count rate to be proportional to the 
energy fiux. This is obviously not strictly correct for a variable 
X-ray slope, possibly leading to inaccuracies in comparing data 
with theory. 
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Figure 7. The F-F correlation in RXTE observations of (a) NGC 
5548, and (b) NGC 4051. Tlie data are from Cliiang & Blaes 
(2003) and Lamer et al. (2003), respectively. The solid lines show 
the best-fit linear dependence, and the dashed lines show the de- 
pendence predicted by eq. j7J with Ep = 180 keV and 140 keV, 
respectively. 



ing a T{F) plot, and it disappears completely in a r(logF) 
space (Fig. 03). 

An important example of broad-band pivoting variabil- 
ity is given by ~ 10^ of 1-day measurements of Cyg X-1 
in the hard state over the 1.5-300 keV energy range (Z02). 
Pivoting is seen here directly in the ASM and BATSE data 
showing a linear r(logF) correlation at low energies and an 
anticorrelation at high energies, not just being inferred from 
a V-F plot. The pivoting variability pattern can also be illus- 
trated by energy-dependent fractional variability. The solid 
and dashed curves in Fig. |U] show the observed fractional 
rms variability as a function of energy is well modelled by 
the theoretical rms due to a variable T with either constant 
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Figure 8. Schematic illustration of power law spectra with 
r = 1.7, 1.9, and 2.1 pivoting at (a) Bp = 10 keV and (b) 
Ep = 200 keV. The two cases correspond to the total X-ray lu- 
minosity being approximately constant and increasing with in- 
creasing r, respectively. The dot-dashed curve shows a possible 
constant reflection component from a distant medium. If the to- 
tal spectrum is fitted by a power law and refiection, the relative 
reflection fraction will increase and decrease with increasing V in 
the cases (a) and (b), respectively. 



Ep and additional energy-independent variability or vari- 
able Ep. We also see that an additional variable component 
is required at the lowest energies, which probably can be 
identified with the observed soft X-ray excess (Ebisawa et 
al. 1996; Frontera et al. 2001a; Di Salvo et al. 2001). Still, 
the 1.5-3 keV Hux remains strongly anticorrelated with the 
100-300 keV flux (Z02). 

Results for Seyfert Is similar to that shown in Fig. |U] 
for Cyg X-1 have been obtained by Markowitz & Edelson 
(2001). They show that the fractional variability in the 2- 
4 keV band in all of 9 Seyfert Is studied by them is higher 
than that in the 7-10 keV band (for long-term measurements 
with 5-day intervals). Given the unknown shape of the de- 
pendence of rms on E, it is not possible to uniquely translate 
their results into the pivot energies. Still, their results are 
compatible with Ep > 10 keV. 

The results above were obtained on time scales from 
days to years for both Seyferts and BH binaries in the hard 
state. Interestingly, the variability patterns in the two cases 
appear qualitatively similar. However, a month time scale in 
Seyferts corresponds to ~ 1 s in BH binaries in units of the 
light travel time accross the gravitational radius. Thus, the 
variability pattern presented here for AGNs physically corre- 
sponds to that on short-time scales in BH binaries. In agree- 
ment with this correspondence, Li, Feng & Chen (1999) find 
the 2-60 keV hardness anticorrelated with the flux (i.e., cor- 
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Figure 9. The rms variability in the hard state of Cyg X- 
1 in one-day averaged measurements by the iJXTiJ/ASM and 
CGRO/BATSE (Z02). The data (crosses) can be modelled by 
either a power-law pivoting of equation IA6i with Ep = 45 
keV, the standard deviation of Ar = 0.23, and an additional E- 
independent variability with era = 0.155 (dashed curve), or a vari- 
able pivot energy (equation lAlSI ) with = 45 keV, Ap = 0.71, 
and Ap = 0.21 (solid curve, almost coincident with the dashed 
one). Note that the non-zero value of the minimum rms of the 
latter model is entirely due to the dispersion of Ep, see Appendix 
A. 



responding to a positive T-F correlation) in the hard state 
of Cyg X-1 at time scales of 10~^ and 1-50 s. Thus, pivoting 
with a high pivot energy is likely to take place in the BH 
binaries on short time scales as well. On the other hand, 
the full picture of variability in the hard state is certainly 
more complicated that this. For example, Li et al. (1999) 
and Wen, Cui & Bradt (2001) find no or weak T-F correla- 
tion in Cyg X-1 on the 0.01-0.1 s and ~ 1-day time scales, 
respectively, whereas it is very strong over hundreds of days 
(Z02). 

Then, spectral state transitions change the character of 
the correlation qualitatively. In particular, the T-F correla- 
tion becomes strongly negative in the soft state of Cyg X-1 
over time scales from 0.01 s to tens of days (Li et al. 1999; 
Wen et al. 2001; Z02). 



4 CORRELATION BETWEEN THE X-RAY 
AND RADIO FLUXES IN BH BINARIES 

A very interesting correlation between radio and X-ray 
fluxes in the hard state of accreting BH binaries has been 
discovered recently (Brocksopp et al. 1999; Corbel et al. 
2000; Gallo et al. 2003; Markoff et al. 2003a). In particu- 
lar, Gallo et al. (2003) show a strong correlation between 
the 15 GHz fiux and the count rate from the ASM from 
daily measurements for 8 BH binaries in the hard state. 
There are two possible origins of this correlation. One is 
that the level of X-ray emission is related to the rate of ejec- 
tion of radio-emitting clouds, forming a compact jet (e.g.. 



Mirabel et al. 1998; Corbel et al. 2000). Another is that 
the X-ray emission of BH binaries is dominated by non- 
thermal emission of the jet (Markoff, Falcke & Fender 2001; 
Vadawale, Rao & Chakrabarti 2001; Markoff et al. 2003a, 
2003b; Georganopoulos, Aharonian & Kirk 2002). 

There are many strong arguments against the second 
interpretation (see also Poutanen & Zdziarski 2003). The 
broad-band X7 spectra of BH binaries in the hard state are 
very well modelled by thermal Comptonization and Comp- 
ton reflection (e.g., Gierlihski et al. 1997; Zdziarski et al. 
1998; Frontera et al. 2001a, b; Wardzihski et al. 2002; Z02), 
see Fig. |5| Apart from those broad-band studies, the ev- 
idence for the presence of Compton reflection from X-rays 
alone is extremely strong (e.g., Gilfanov et al. 1999; CGROO; 
Revnivtsev et al. 2001; Nowak et al. 2002). This implies that 
the X-ray emission is not strongly beamed away from the 
disc. 

The thermal- Compton origin of the primary X-ray emis- 
sion is strongly supported by a remarkable uniformity of 
both the energy (within a factor of ~ 2) and shape of the 
high-energy cutoffs of BH binaries in the hard state ob- 
served by OSSE (Grove et al. 1998). This cutoff is natu- 
rally accounted for by thermostatic properties of thermal 
Comptonization as well as pair production (e.g., Malzac, 
Beloborodov & Poutanen 2001), as it corresponds to the 
transition to relativistic temperatures. At higher tempera- 
tures, cooling becomes extremely efficient and copious pair 
production starts. This reduces the energy available per par- 
ticle causing the temperature to decrease. 

In the framework of the model with synchrotron emis- 
sion of non-thermal electrons accelerated in the first-order 
Fermi process in a shock, the maximum energy of the syn- 
chrotron photons is given by balance of the acceleration and 
energy loss time scales. This gives the critical synchrotron 
energy corresponding to the maximum Lorentz factor of the 
accelerated particles, 7max, as 



2 o4 ■ r,2 

nieC S sma p^h 



Qf 



25 



(8) 



where /3shC is the shock velocity, ^ (satisfying 1 //3sh ^ C 1) 
gives the efficiency of the shock acceleration, at is the fine- 
structure constant, a is the pitch angle (with (sin a) ~ 1/2) 
and nicC^ is the electron rest energy. We see that for f3^i-^/^ ~ 
10""^ we can reproduce the high energy cutoff of BH binaries 
at ~ 100 keV, as pointed out by Markoff et al. (2001). 

However, even with the fine-tuning of /Sfh/Ci the cutoff 
in the electron distribution is unlikely to be sharp. One ef- 
fect is that of electron energy loss, which makes the steady- 
state distributions to cut off gradually (e.g.. Kirk, Rieger 
& Mastichiadis 1998). Also, the conditions in the acceler- 
ating region (determining the efficiency parameter, ^) are 
unlikely to be completely uniform. Under many astrophys- 
ical circumst ancles, the shape of the tail of a distribution 
is given by an exponential cutoff (e.g., in a Maxwellian). If 
we assume the power law electrons, N{'y) = 1)7"*', are cut 
off as exp(— 7/7inax), the resulting synchrotron spectrum, 
Fe = EdN/dE, is given by. 



Fe 



K 



i+p 

dyy 2 E2+P 



K5/3(y), (9) 



where e — EjE^ 



E„ is the exponential integral of the 
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order n, K„ is the modified Bessel function of the second 
kind, the normahzation constant is, 



K = 



sin a 



(10) 



B is the magnetic field strength, e is the electron charge, h 
is the Planck constant, and Eq is the constant factor in the 
critical synchrotron energy, Ec, corresponding to emission 
by electrons with the Lorentz factor 7, 

^ 2 3ehB sin a ,,,, 
Ec = , Eo = — . 11 

The spectrum of equation Q exhibits a rather gradual cut- 
off, in fact much more gradual than that of a typical hard- 
state spectrum of Cyg X-1, as as illustrated in Fig.|2| Indeed, 
the synchrotron model of the hard state of Cyg X-1 shown 
in fig. 3a in Markofi' et al. (2003b) overestimates the 1 MeV 
flux (McConnell et al. 2002) by a factor of 8 when matched 
to its 100 keV flux. 

We have also studied various other forms of the elec- 
tron high-energy cutoff. We have found that anything more 
gradual than a sharp cutoff at 7max, when 



f^)K5/3(y), (12) 



gives photon spectra with cutoffs not strong enough to 
match the cutoff of Cyg X-1, as illustrated in Fig.|21 

A related issue is the spectral index of the spectrum. 
The model of Markoff et al. (2001) relies on particle accel- 
eration in shocks. As they note, the accelerated electrons 
have then a power law index > 1.5, where 1.5 corresponds 
to highly relativistic shocks. As the terminal bulk motion of 
jets in BH binaries is only mildly relativistic, with the typ- 
ical bulk Lorentz factor of a few, it is unlikely that a shock 
close to the base of a jet is highly relativistic. Electron en- 
ergy losses steepen the steady-state electron distribution by 
unity (e.g.. Kirk et al. 1998), to an index p > 2.5, and the 
corresponding photon index is F = (p -|- l)/2. Markoff et al. 
(2001) have been able to fit this model to the spectrum of 
XTE J1118-I-480, which has F ~ 1.8. However, many BH bi- 
naries have X-ray spectra harder than that (see Fig.0. The 
resulting discrepancy is illustrated in Fig. |21 where we as- 
sumed the limiting p = 2.5 in the synchrotron model spectra 
shown. 

In order to deal with this rather serious problem, 
Markoff et al. (2003a) assumed that electron acceleration is 
continuous through the jet rather than confined to a shock. 
Then, the steady-state electron index equals that of the ac- 
celeration process (e.g.. Kirk et al. 1998). Markoff et al. 
(2003a) do not specify details of this acceleration process. A 
natural candidate seems to be here the second-order Fermi 
acceleration in a turbulent medium (e.g., Jones 1994), in 
which case the electron index is constrained only by p > 1. 
This scenario is possible in principle, but then the accelera- 
tion depends on the unknown properties of the accelerating 
medium and equation ||HJ , determining the high-energy cut- 
off, no longer applies. 

The problems with the energy and shape of the high- 
energy cutoff are even more severe for the non-thermal 
Compton model (Georganopoulos et al. 2002). It requires 
that the product of the seed photon energy and the square of 
the maximum electron energy is fine-tuned, the high-energy 



cutoff in the electron distribution very sharp, and the seed 
photon distribution is very narrow. The last condition is not 
fulfilled in any of the models of Georganopoulos et al. (2002) , 
which rely on scattering of either disk or stellar blackbody 
photons. Furthermore, the X-ray jet size in the model of 
Cyg X-1 of Georganopoulos et al. (2002) is 10^^ cm, which 
is clearly inconsistent with the peak of the power spectrum 
per logarithm of frequency being at some occasions at ~ 10 
Hz (e.g., Revnivtsev, Gilfanov & Churazov 2000). 

As discussed in Section |5| the amplitude of Compton 
refiection and the Fe Ka flux (Section |SJ imply that dense 
and rather cold material occupies a solid angle of f2 ~ tt as 
viewed from the X-ray source. Smearing of these components 
(e.g., Frontera et al. 2001a; Di Salvo et al. 2001; GCROO) 
and their correlation with F (Sections |5| [SJ clearly identify 
the reflector with the accretion disc (Section]^ and imply 
most of the X-ray source is within ~ 30-100 gravitational 
radii from the BH. On the other hand, models of Markoff et 
al. (2001, 2003a, 2003b) and Georganopoulos et al. (2002) 
ignore that spectral components. 

Yet another piece of evidence against a substantial part 
of X-rays being non-thermal is provided by spectral vari- 
abflity. In the case of Cyg X-1, the ASM/BATSE data show 
a pivot around ~ 40-50 keV (Z02). This is consistent with 
both the X-ray hardness-flux data and the rms variability 
strongly increasing with decreasing X-ray energy, see Ta- 
ble 1 and Fig. |5] respectively. The characteristic Ar from 
those data is ~ 0.2-0.3. This power-law spectral variabil- 
ity when extended to the turnover energy (when the syn- 
chrotron source becomes optically thick) at ~ 1 eV (e.g., 
Markoff et al. 2001, 2003a) would imply variability below 
this energy by a factor of ~ 10^. However, the range of the 
variability of the 15 GHz flux correlated with the ASM flux 
in Cyg X-1 is only by a factor of several, basically the same 
as the range of the variability of the ASM flux itself (Gallo 
et al. 2003), contradicting the synchrotron origin of X-rays. 

If both radio and X-rays were indeed due to non-thermal 
synchrotron emission (Markoff et al. 2001, 2003a, 2003b), 
their observed variability pattern should yield the X-ray rms 
virtually independent of energy. Then, in the model with 
non-thermal Comptonization of photons from the compan- 
ion star (Georganopoulos et al. 2002), the seed photon flux is 
just constant, and the fractional variability should increase 
with photon energy globally. These model predictions are in 
strong disagreement with the data shown in Fig. |^ and, in 
particular, with the ASM 1.5-3 keV flux being stronly an- 
ticorrelated with the 100-300 keV flux from BATSE (Z02). 
A similar anticorrelation occurs in Cyg X-3 (McCoUough et 
al. 1999b), and, in fact, the BATSE flux in that object is 
anticorrelated with radio flux (McCoUough et al. 1999a). 

On the other hand, if we extrapolate the rms depen- 
dence of Fig.|5]to low energies and make a plausible assump- 
tion of the fractional rms being < 1, we obtain the charac- 
teristic energy of a fraction of keV. This is in very good 
agreement with the temperature of seed photons for ther- 
mal Comptonization observed in Cyg X-1 to be kThh ~ 0.15 
keV (Ebisawa et al. 1996; Di Salvo et al. 2001). 

All these results strongly support the interpretation of 
the correlated radio emission as being due to ejection of 
clouds from the X-ray source, which could be similar to coro- 
nal mass ejections observed at the Sun. This interpretation 
is also supported by observed time lags of the radio emission 



© 2003 RAS, MNRAS OOO.IIHT^ 



12 A. A. Zdziarski et al. 



with respect to the X-rays, e.g., in GRS 1915+105 (Mirabel 
et al. 1998). In fact, the base of the jet may be just the hot 
inner flow (e.g. Fender 2002). However, the electron distri- 
bution in that base is still thermal, as argued above. 

Finally, we note a problem with internal consistency 
in the X-ray jet model of GRS 1915-1-105 of Vadawale et 
al. (2001). Namely, integrating the model non-thermal syn- 
chrotron emission shown in their fig. 4 yields the total jet 
luminosity of ~ 10*^ erg s~^ (at a distance of 12 kpc, jet 
velocity of 0.9c, and i — 70° adopted in that paper; note 
that this emission is beamed away from the observer). On 
the other hand, the total jet power given in that paper is 
4.3 X 10=*^ erg s"\ This power includes the proton rest mass, 
and at the velocity of 0.9c, the kinetic power is ~ 2 x 10^^ 
erg s^^. Thus, the jet radiative luminosity is ~ 50 times the 
kinetic power, which strongly violates energy conservation. 
Also, the radiative luminosity is usually much less than the 
kinetic power (unless the jet stops entirely, in which case 
only the two luminosities can be comparable), which fur- 
thermore increases the energy conservation problem. 



5 THEORETICAL INTERPRETATION 

5.1 Feedback between cold and hot media 

The F-F and il-F correlations occur likely due to interac- 
tion between cold and hot media. The former correlation 
(also manifesting itself as pivoting) is likely due to variabil- 
ity in the flux/luminosity of seed soft photons (emitted by 
some cold medium) irradiating a hot plasma being much 
stronger than the variability of the flux/luminosity from the 
hot plasma itself. The main radiative process in the hot 
plasma needs to be thermal Compton upscattering of the 
soft photons. Then, the larger the irradiating flux of seed 
photons, the softer and stronger the X-ray spectrum. This 
variability pattern in the case of a constant hot-plasma lumi- 
nosity is shown, e.g., in fig. 3 of ZGOl and fig. 14 of Z02. The 
pivot point is then somewhere in the middle of the broad- 
band spectrum. If the hot-plasma luminosity increases as 
well but slower than the irradiating flux, the pivot is at high 
energies. 

The emission of the cold medium irradiating the hot 
plasma may be partly (or wholly) due to reprocessing of 
the emission of the hot plasma. Then, the increased cool- 
ing of the hot plasma (which softens the X-ray spectrum) 
is associated with more Compton reflection (accompanying 
reprocessing by a Thomson-thick medium). This gives rise 
to an fi-F correlation (ZLS99; CGROO). 

The two patterns can occur together or independently. 
If the variability of soft seed photons is intrinsic (not from 
reprocessing), a, T-F correlation will not be accompanied 
by an Q-T one. This may happen in narrow-line Seyfert Is. 
On the other hand, the variable seed photons may be due 
to reprocessing but from emission of a highly variable hot 
plasma. Then, there will be an fi-F correlation but not a F- 
F one, which is the case for, e.g., Cyg X-1 on intermediate 
time scales (Gilfanov et al., in preparation). 

A number of specific geometries have been proposed. 
In one, there is a variable radial overlap between the hot 
and cold accretion discs (ZLS99; Poutanen, Krolik & Ryde 
1997). In another, the cold disc extends all the way to the 



minimum stable orbit and the hot plasma forms a corona 
with a mildly relativistic velocity directed either away from 
the disc or towards it (Beloborodov 1999a, 2001; Malzac et 
al. 2001). Also, static coronae have been considered in the 
context of spectral correlations (e.g., Haardt, Maraschi & 
Ghisellini 1997; MFOl). 

5.1.1 Variable overlap between hot and cold flows 

ZLS99 have interpreted the il-F correlation as due to feed- 
back in an inner hot (thermal) accretion flow surrounded by 
an overlapping cold disc. Then, the closer to the central BH 
the cold disc extends, the more cooling of the hot plasma 
by blackbody photons (both reprocessed and from intrinsic 
dissipation), and the softer the spectrum. The effect of the 
cooling by the UV photons of the X-ray emitting plasma is 
seen, e.g., in NGC 7469, where there is a positive correlation 
between the UV flux and F (NOO). Also, the F-F correlation 
(Section El indicates the dominant effect of plasma cooling 
by a variable seed-photon flux. At the same time, the cold 
disc subtends a larger solid angle from the point of view 
of the hot plasma, and thus there are more reflection pho- 
tons in the spectrum. The solid curve in Fig. shows the 
prediction of a simple version of this model (ZLS99). 

This model also naturally accounts for the correlations 
of the reflection strength with both the QPO frequency and 
the degree of relativistic smearing seen in BH binaries (Gil- 
fanov et al. 1999; GCROO; Revnivtsev et al. 2001). The cor- 
relation with the QPO frequency is expected because that 
frequency is very likely to be related in some way to the Ke- 
plerian frequency at the inner edge of the cold disc, which 
increases with the decreasing disc inner radius. At the same 
time, the closer the reflecting medium to the BH the higher 
the degree of the relativistic smearing. 

The model naturally explains the pivoting behaviour of 
the X-ray spectra (manifesting itself in X-rays as a F-F cor- 
relation) as driven by the variable flux of irradiating seed 
photons (ZGOl; Z02). Also, the speciflc F-F correlations ob- 
served in many accreting BHs indicate the bolometric lu- 
minosity increases in those sources with the increasing F 
(Section OJ. This scenario can account for this behaviour if 
the inner radius of the cold disk decreases with the increas- 
ing accretion rate. Such a behaviour is indeed postulated in 
models of advection-dominated accretion (e.g., Esin et al. 
1997) as well as predicted by models of accretion disk evap- 
oration (Meyer, Liu & Meyer-Hofmeister 2000; Rozahska & 
Czerny 2000). On the other hand, it also appears likely that 
some instabilities can affect the disc truncation radius even if 
the bolometric luminosity remains approximately constant, 
which is observed in some objects (see Section 

Chiang & Blaes (2001, 2003) and Chiang (2002) have 
shown that detailed versions of this model can also explain 
the overall optical/UV/X-ray variability in a few Seyferts 
(NGC 3516, NGC 5548, NGC 7469). Note that their calcu- 
lations require the variable overlap to be mostly achieved by 
the radius of the hot plasma being variable. 

We note that this model does not readily explain 
n/27r > 1 sometimes observed (see Figs. 13 |lji). Further- 
more, scattering of the reflected photons in the hot flow 
will further reduce the observed fl (see a discussion in Be- 
loborodov 2001). These problems can be possibly resolved 
by anisotropy of the emission of the hot plasma or the outer 
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disc being concave. It is also possible that detections of such 
large reflection are due to imperfection of the spectral mod- 
els used and/or data inaccuracies. 

In the framework of this model, radio emission arises 
from outflows in the hot inner flow (e.g., Blandford & Begel- 
man 1999). Such a scenario is described, e.g., in Fender 
(2002). Note that the jet emission in BH binaries is usually 
quenched in the soft state (Brocksopp et al. 1999; Corbel et 
al. 2000) , in which the hot inner flow most likely disappears 
and is replaced by a hot corona (e.g., Gierlihski et al. 1999). 

5.1.2 Static coronae 

Haardt et al. (1997) have studied spectral correlations in 
a static disc corona model. They found that if the corona 
is dominated by pairs, the 2-10 keV spectral index, F 
is rather insensitive to changes of the flux in the same en- 
ergy range, with typical AF ~ 0.2 for a change of the flux 
by 10. This is clearly much less than the observed spectral 
variability in many sources (Section]^, which rules out this 
model. The predictions for coronae not dominated by pairs 
depend, in turn, on the choice of the coronal optical depth, 
r. At r>0.3, F would decrease with the flux, contrary to 
the data, but coronae with lower r could be reconciled with 
data. 

Then, MFOl found that static patchy coronae could re- 
produce the observed F-_F correlations if the luminosity of 
an active region increases with its increasing size at a given 
height. The increased size increases the feedback of soft ra- 
diation from the disc, which makes, in turn, the spectrum 
softer. MFOl and Shih et al. (2002) have shown that this 
model fits well the F-F data for MCG -6-30-15 from RXTE 
and ASCA, respectively. 

On the other hand, this model yields an Q-T anticorre- 
lation rather than correlation because an increased size at 
a given height (leading to the softening of the spectrum, 
see above) increases the degree of obscuration of the re- 
flected radiation (Beloborodov 1999b, Malzac et al. 2001). 
This property is also shared by static disc coronae in general, 
and it follows from the dependencies shown by Haardt et al. 
(1997). Also, static coronae cannot by themselves explain 
the formation of radio jets. 

5.1.3 Dynamic coronae 

An interpretation of the Sl-F correlation alternative to the 
variable overlap of two accretion flows (Section I5.1.H is 
that with mildly relativistic coronal inflows/outflows (Be- 
loborodov 1999a; Malzac et al. 2001). The higher the speed 
of the coronal outflow, the less feedback with the underly- 
ing disc and the harder the spectrum. Inflows can, in turn, 
account for i}/2n > 1. 

It is not clear how to account for the F-F correlation in 
this model. This correlation requires that the observed lu- 
minosity of the hot plasma either stays constant or increases 
with a decrease of the outflow velocity. On the other hand, 
an opposite behaviour of constant luminosity of the seed 
photons and variable luminosity of the hot plasma (which 
results in a T-F anticorrelation) is seen in the soft state of 
Cyg X-1 (Churazov, Gilfanov & Revnivtsev 2001; Z02), in 
which case coronal models are widely accepted (e.g., Pouta- 
nen et al. 1997; Gierlihski et al. 1999; Churazov et al. 2001). 



In the framework of this model, radio emission may arise 
from the jet being formed by the coronal outflow (Merloni & 
Fabian 2002). However, the quenching of the radio emission 
in the soft state (Brocksopp et al. 1999; Corbel et al. 2000) 
is not readily explained in this model. 



5.1.4 Further physical implications 

The feedback models explain the fl-V correlation in terms 
of Comptonization of blackbody photons by a hot thermal 
plasma. The effective reflection solid angle in this model 
can be linked to the geometry of the source, which will also 
determine the amplification factor. A, of the Comptoniza- 
tion process. On the other hand, the spectral index follows 
from the energy (and pair) balance. Since the character- 
istic blackbody seed photon energy is much lower in AGNs 
than in BH binaries, a generic prediction of these models is 
that for a given Q, the value of F will be higher in AGNs 
than that in BH binaries (Beloborodov 1999b; Malzac et al. 
2001). This is illustrated in Fig. llOl fwhich calculations were 
performed using a code of Coppi 1999). When the plasma 
compactness (oc L/R, where L and R are the luminosity 
and size, respectively), is low, pair production is negli- 
gible and the energy balance is satisfied by adjusting the 
electron temperature, kT, see Fig. llOb . On the other hand, 
at high compactness, the energy balance is achieved by ad- 
justing the Thomson optical depth of the pairs, see Fig. llOb . 
(In the calculations shown in Fig. llOb . we assumed the pres- 
ence of a weak high-energy power-law tail, containing 0.05 
of the total energy beyond the Maxwellian electron dis- 
tribution, as suggested by the results of McConnell et al. 
2002.) The calculations shown in Fig. 1101 with F changing 
from 1.63 to 1.85 at the assumed f2/27r = 0.4, explain well 
the difference between the T-Q correlations for the BH bi- 
naries (Fig.EJl and AGNs (Fig. El . 

An important issue here is also the possible presence of 
a correlation between F and kT for a given feTbb- The feed- 
back models require that the parameters of the hot thermal 
plasma adjust to variable cooling to satisfy the energy bal- 
ance. If T is constant, this can be done only by changing 
kT. Then, an increased cooling yielding an increase of F 
is associated with a decrease of kT. On the other hand, it 
is possible that t also changes, e.g., responding to changing 
accretion rate, which would produce a more complicated be- 
haviour (e.g., Chiang & Blaes 2001, 2003). Also, if pairs 
are present, the main adjusting parameter is r. In fact, in 
the example in Fig. llOb . the softening of the spectrum is 
associated with an increase (rather than decrease) of kT. 

So far there are only limited observational data on cor- 
relation between F and the high-energy cutoff. P02b found 
a positive correlation between F and the e-folding energy in 
their sample of Seyfert Is, which would be contrary to expec- 
tations of the simplest model with constant r and the vari- 
ability of F due to changing kT. However, MP02 found that 
this correlation may be an artifact of the e-folded power- 
law model when used to fit thermal-Compton spectra. They 
simulated actual thermal Comptonization spectra at a con- 
stant kT = 100 keV and a range of r, and then found that 
the simulated spectra when fitted by an e-folded power law 
yield a correlation resembling that of P02b. 
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Figure 10. Illustration of the effect of changing fcTbb of the seed 
photons from 5 eV (dashed curves; characteristic to AGNs) to 
150 eV (solid curves, characteristic to BH binaries). In both cases 
the Comptonization amplification factor, A = 10. The plasma is 
dominated by (a) electrons and protons (e~p) and (b) pairs. 
The e~p plasma parameters are t = 1.5, kT = 57 keV, and r = 
1.5, kT = 87 keV, in the AGN and BH binary cases, respectively, 
which yields the corresponding 2-10 keV spectral indices of F = 
1.85 and 1.63. The plasma parameters are r = 1.0, kT = 84 
keV, and t = 1.6, kT = 72 keV, in the AGN and BH binary cases, 
respectively. In all cases, Q/2it = 0.4 at i = 30°. 

5.2 Reflection from strongly ionized disc 

Done & Nayakshin (2001) have shown that the optical depth 
of a highly ionized surface layer on top of a strongly irra- 
diated accretion disc increases with the increasing hardness 
of the irradiating spectrum. This gives rise to an appar- 
ent fl-F correlation as there is less unscattered reflection for 
harder spectra (although the actual reflection solid angle is 
Q = 2tt = constant). A separate mechanism is then needed 
to account for the T-F correlation, similarly to the case of 
the coronal outflow model. 

A diagnostic that could yield the actual solid angle of 
the reflector regardless of the ionization level is a measure- 
ment of reflection at high energies. Namely, the reflection 
spectrum is cut off at high energies with a functional form 
independent of the ionization state (White, Lightman & 
Zdziarski 1988). Fig. 1111 illustrates this point for an inci- 
dent spectrum from Comptonization (Poutanen & Svensson 
1996) with kT = 120 keV and r = 2 in spherical geometry, 
yielding the 2-10 keV index of F = 1.5. This F corresponds 
to the case with the lowest fitted neutral reflection from a 
strongly ionized medium with the actual f2/27r = 1 in Done 
& Nayakshin (2001). To illustrate this effect, we assumed 
that the reflecting medium is so strongly ionized that vir- 
tually no Fe K edge appears in the reflection component 
(dotted curve in Fig. IIH . Although the reflected component 
is simply a power law at low energies, it still does have a high 
energy cutoff due to Klein-Nishina effects. The solid curve 
shows the total spectrum. Then, the dashed curve shows the 
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Figure 11. Comptonization spectrum with negligible reflection 
(dashed curve) compared to a spectrum with f!/27r = 1 and the 
reflecting medium being almost fully ionized (solid curve). The 
dotted curve shows the ionized reflection component alone. We 
see that the form of the high-energy cutoff in the reflection spec- 
tral component leads to the total spectrum being distinctly differ- 
ent from the Comptonization spectrum alone (although the two 
coincide at low energies). See Section 15.21 for details. 

incident spectrum without reflection normalized to coincide 
with the spectrum from the strongly ionized medium. We 
see that although the spectra at < 20 keV are indeed barely 
distinguishable, the form of their high-energy cutoffs is very 
different. Thus, this model can be tested using broad-band 
X7 spectra extending to several hundred of keV. 

A set of ^ 1-1000 keV spectra of Cyg X-1 in the hard 
state was analyzed by Gierlihski et al. (1997), who found 
the form of the high-energy cutoff compatible with thermal 
Comptonization. Although those authors have not consider 
the ionized-reflection model, those spectra have not shown 
any hint of the softening of the cutoff of the type seen in the 
solid curve in Fig. 1111 rather they indicated some additional 
hardening (interpreted by them as due to the presence of an 
additional spectral component from saturated Comptoniza- 
tion). 

Independently of the present work. Barrio, Done & 
Nayakshin (2003) have recently fltted the ionized reflec- 
tion model to broad-band spectral data from the PCA and 
HEXTE detectors of RXTE for Cyg X-1. They have found 
this model to be ruled out by the lack of a break in the 
HEXTE data, in favour of the model with a truncated disk 
and intrinsically weak reflection. Future tests of this model 
are highly desirable, e.g., using data from INTEGRAL and 
ASTR0-E2. 

The effect shown in Fig. Illl also demonstrates that the 
albedo, a, of a fully ionized medium energy-integrated over 
a spectrum characteristic to accreting BHs is still substan- 
tially less than unity (see also White et al. 1988). Indeed, 
in the specific case in Fig. 1111 a = 0.48, only moderately 
larger than the corresponding albedo for a neutral refiector, 
a — 0.34. This shows that a fully ionized surface layer of 
an accretion disc would radiate substantial soft thermal fiux 
even in the absence of an internal dissipation. This may rule 
out the ionized-disc model of the hard state of Cyg X-1 of 
Young et al. (2001). In general, the resulting constraints on 
the geometry will be relatively similar to those in the case of 
a nearly-neutral disc (e.g.. Stern et al. 1995), and, in partic- 
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ular, will still rule out the presence of a homogeneous corona 
above the disc. 

5.3 Time lags and reflection 

NOO have pointed out that the presence of a constant re- 
flection component from a distant medium (presumably a 
molecular torus) in AGNs together with a variable power 
law (from hot plasmas in the vicinity of the BH) can lead to 
an apparent correlation between F and Q if reflection in the 
fitted model is tied to the variable power law. In particular, 
when the power law varies with a characteristic pivot energy 
of < 10 keV, the fitted il-F dependence will be positive (i.e., 
the softer the spectrum, the larger the relative contribution 
of reflection), as illustrated in Fig.|S^. 

This idea has been developed in detail by MP02. They 
have considered cases with the pivot energy of 2, 5 and 10 
keV, and have shown that the fl-T dependence in Seyferts 
(Fig-EJ can be qualitatively explained by this type of model. 
However, we have shown in Section |21 that the pivot energy 
in many Seyferts is > 100 keV. Then, if this effect were 
important for a large fraction of Seyferts, an Q-T anticorre- 
lation (see Fig.lHJj) would be observed for a large fraction of 
sources, in conflict with the data. This rules out the origin 
of the global fi-F dependence in Seyferts from this effect. 
On the other hand, this effect may be important in some 
Seyferts, and may explain Q-T anticorrelation in repeated 
observations of certain objects (Done et al. 2000; P02a). 

As noted by MP02, this effect cannot play any role in 
BIf binaries. Thus, this model requires completely different 
origin for the Q-T correlations in those objects and in AGNs. 



6 THE Fe Ka EMISSION 

The results presented in Section |21 were mostly obtained by 
fits to the Compton-reflection continuum. Ifowever, the pro- 
cess of bound-free absorption in the reflecting medium is 
usually followed by emission of a fluorescent Fe Ka line, 
which equivalent width, Wko , can be tied to that of the re- 
flection (George & Fabian 1991; Zycki & Czerny 1994). Con- 
versely, if relativistic broadening is seen in the line emitted 
by the reflecting medium, the same broadening should be 
seen in the reflection spectrum. Thus, the measured param- 
eters of the Fe Ka line should be in agreement with those of 
reflection. Such an agreement appears to be seen in BH bina- 
ries (Gilfanov et al. 1999; GCROO; Revnivtsev et al. 2001). 
Some disagreement in GX 339-4 claimed by Nowak et al. 
(2002) appears to be largely caused by a diffuse line com- 
ponent from the interstellar medium dominating in X-ray 
weak states of this source (Wardzihski et al. 2002). 

On the other hand, the situation is more complex in 
Seyferts. Unlike BH binaries, Seyferts possess signiflcant Fe 
Ka line component from distant matter, e.g., the broad- 
line regions or molecular torii. Often, that distant matter is 
Thomson thin, and then it gives rise to the line but not to an 
Fe K edge, see, e.g., section 7.1 and fig. 13 in Wozniak et al. 
(1998). As shown in that figure, the optically-thin reflection 
component is then relatively negligible. Then, the presence 
of two line components, one constant (with or without asso- 
ciated reflection) and one variable in addition to a variable 
intrinsic component may cause a rather complex behaviour 



of variability. For example, the presence of constant line and 
reflection and a power law pivoting at ~ 10 keV leads to an 
anticorrelation of Wkoi and ^1/2tt measured with respect to 
the power law (see Fig.|H^, also MP02). Such an effect might 
have been seen by Chiang et al. (2000) in NGC 5548. 

Given this complexity, it is important to self- 
consistently account for various components of the line and 
reflection, their relativistic broadening, and their connec- 
tion to the incident spectral component. A self-consistent 
treatment should in general include at least two line compo- 
nents, narrow and broad, with the broadening (e.g., Gaus- 
sian or disc-like) and redshift of the latter applied in the 
same way to the reflection (e.g, Zycki et al. 1998, 1999; 
GCROO, Lubiriski & Zdziarski 2001; Zdziarski et al. 2002a). 
The strength of reflection should then be tied to that of 
the line (George & Fabian 1991; Zycki & Czerny 1994). The 
dramatic effect of inconsistent treatment of the line and re- 
flection is shown for the case of the disc broadening (Fabian 
et al. 1989) in fig. 3 of Zdziarski et al. (2002a). 

Then, Shih et al. (2002) showed that Fe Ka line fiux in 
MCG -6-30-15 remained approximately constant with vary- 
ing 3-10 keV ASCA count rate. However, the formation of 
the line is due to photons only above 7.1 keV or more. In the 
case of MCG -6-30-15, F ~ 2, and then the fraction of the 
measured ASCA counts above 7 keV in their data is only 
~ 0.03. The fraction of counts above 8 keV, appropriate for 
a moderately ionized disc, possibly present in MCG -6-30-15 
(e.g., Ballantyne & Fabian 2001), is as low as ~ 0.01. Thus, 
the lack of a correlation between the 3-10 keV counts and 
the line photon fiux cannot be taken as conclusive evidence 
against the Fe Ka photons being due to irradiation by pho- 
tons from the observed continuum. Similarly, the fraction 
of counts above 7 and 8 keV in the MCG -6-30-15 data of 
VEOl is only ~ 0.2 and ~ 0.1, respectively. The data of 
VEOl imply a pivot energy at about 80 keV (Table 1); thus, 
the variability of the fiux above ~ 7-8 keV will be in gen- 
eral less than that in the 2-10 keV band used by them (as 
also found by Markowitz & Edelson 2001 for Seyfert Is in 
general) . 

As noted in Section |21 (as well as by MP02) , there are 
many AGNs with the pivot at < 10 keV (e.g., 3C 120). Then, 
the 7-10 keV photon flux may be anticorrelated with the 2- 
10 keV one. Even if the pivot energy is in general at ^ 10 
keV, occasionally the same object may show a pivot at < 10 
keV, as in the case of NGC 5548 (Nicastro et al. 2000). 
Thus, it is important to use the band above the Fe K edge 
in studies of the continuum-line correlations. 

Self-consistent treatments of the line and reflection was 
applied, e.g., to the Seyfert galaxy IC 4329 A by Done et 
al. (2000), who found no disagreement between the line and 
reflection. A similar method was applied to NGC 4151 by 
Zdziarski et al. (2002a) , who also found the broad-line com- 
ponent compatible with reflection, in spite of previous claims 
to the contrary (e.g., Wang, Zhou & Wang 1999). Then, 
Lubinski & Zdziarski (2001) have applied the same method 
to average ASCA spectra of Seyfert Is and found the depen- 
dence of the Wkq of the broad component of Fe Ka lines 
to obey the dependence on F compatible with that of the 
reflection continuum. Although Yaqoob et al. (2002) have 
pointed out that the statistical weights used by Lubinski 
& Zdziarski (2001) were not justifled statistically, the effect 
of their correction appears to be only minor (Lubinski, in 
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preparation). Indeed, the results of P02b for Seyfert Is ob- 
served by BeppoSAX show an increase of Wkci with F very 
similar to that obtained by Lubihski & Zdziarski (2001). 



7 CONCLUSIONS 

We have considered correlations between various spectral 
properties of accreting BHs in Seyfert galaxies and X-ray 
binaries, with particular emphasis on the correlations be- 
tween the X-ray spectral index, strength of Compton reflec- 
tion and the X-ray flux. The main results can be summarized 
as follows. 

Using published data of the observations of Seyferts 
with Gmga, RXTE and BeppoSAX, we have critically re- 
evaluated the evidence for presence of correlation between 
the X-ray spectral index and strength of Compton reflection. 
We conclude that when considering a large number of obser- 
vations of a large sample of objects, the existence of global 
correlation between these two parameters is established be- 
yond any reasonable doubts. 

Smallness of the error bars in comparison with the ex- 
tent of the correlation and good agreement of the results 
obtained by various satellites confirm that the correlation 
cannot be an artifact caused by statistical or systematic ef- 
fects. The ratios of the spectra with different values of reflec- 
tion in BH binaries demonstrate that the correlation cannot 
be a consequence of a trivially inadequate spectral model. 
We note, however, that the particular values of the spectral 
index and, especially, of the strength of the Compton reflec- 
tion, do depend on the details of the spectral approximation. 
This fact should taken into consideration when comparing 
results obtained by difi'erent authors. 

The Q.-V correlation shows significant spread, larger 
than the statistical uncertainties of the data. It is not clear 
at present to which degree this spread is intrinsic to the 
sources and to which degree it is due to imperfectness of 
the spectral model and/or difi'erence in the details of the 
spectral approximation. 

Distinction should be made between classes of objects 
and multiple observations of individual objects. In the case 
of luminous BH binaries, spectral variability of individual 
sources obeys the general Q.-T correlation obtained for BH 
binaries as class. In the case of Seyfert galaxies, the cor- 
relations holds for a class of objects but can be violated 
for repeated observations of individual objects, remaining, 
however, within the spread of the global correlation. In our 
opinion, these results are still inconclusive and require fur- 
ther investigation as a number of complications are involved 
in the case of Seyfert galaxies in comparison with X-ray bi- 
naries. The most obvious among those are lower statistics 
due to significantly lower brightness of Seyferts and exis- 
tence of molecular tori and broad line regions, which can 
give additional contribution to the reflected component, and 
uncorrelated with the Comptonized emission on short time 
scales. 

The physical interpretion of the Q.-T and T-F correla- 
tions will advance our understanding of the geometry of the 
accretion flow in X-ray binaries and Seyfert galaxies and 
impose valuable constrains on the theoretical models. At 
present, the correlations appear to be a natural consequence 
of co-existence of cold media (e.g., an accretion disc) and a 



hot Comptonizing cloud in the vicinity of the BH. Their ge- 
ometrical closeness results in double feedback between the 
two components of the accretion flow. The cold medium pro- 
vides seed soft photons for the Comptonization as well as 
it reprocesses and reflects the hard radiation from the hot 
plasma. If significant fraction of soft seed photons is due 
to reprocessing of the hard radiation of the hot cloud, the 
more reprocessing results in more cooling of the hot plasma 
and, correspondingly, in the softer X-ray spectra. Among 
the several specific geometries proposed, the most promising 
appears to be the disc-spheroid model with variable overlap 
between the hot and cold components of the accretion flow. 

Since the characteristic temperature of the accretion 
disk in AGNs is much lower than in BH binaries, a generic 
prediction of this type of models is that for a given value 
of reflection the value of the spectral index will be higher 
in AGNs than in BH binaries. This prediction is in good 
agreement with the observed behavior. 

We presented a diagnostic to test an interpretation 
of the Q,-T correlation as due to strong ionization of the 
disc surface layer. It utilizes the independence of the Klein- 
Nishina cutoff in the reflected spectrum of the ionization 
state. Applied to existing data, this diagnostic does not sup- 
port that interpretation. 

We found that the pattern of broad-band spectral vari- 
ability of Seyfert galaxies on day-to-month time scale in- 
cludes a pivoting of a power law spectrum with the pivot 
at a high energy, usualUy above a few hundred keV. The 
pivoting well explains the linear correlations between the 
logarithm of the X-ray flux and F seen in both BH binaries 
and Seyferts. Then, the observed high pivot energies rule 
out the interpretation of the Q,-V correlation is Seyferts as 
due to a time-lag effect. 

We discuss also the correlation between X-ray and ra- 
dio fluxes in BH binaries and its physical implications. We 
conclude that this correlation is most likely due to relation 
between the level of X-ray emission and the rate of ejection 
of radio-emitting clouds forming a compact jet. Although 
peculiar sources might exist, it seems highly unlikely that 
the correlation is due to synchrotron origin of the X-ray 
emission from BH binaries in general. 
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APPENDIX A: PIVOTING 

Let us consider spectral variability consisting of pivoting. 
Namely, an initial spectrum (a power law or not) is multi- 



plied by {E/Ep) 



/, where Ep is the pivot energy and 



5 is a perturbation spectral index. In particular, when the 
initial spectrum is a power law, we have a variable power- 
law photon spectrum, AN /AE = C(E / Ep)^^ , where C is 



a constant. Then there is linear relation between the loga- 
rithm of the monochromatic flux at a given energy, E, and 
the variable slope. 



ln{AN/AE) = ln(£;p/£;)r + InC. 



(Al) 



A similar linear relation holds approximately for the energy 
flux in an interval from Ei to E2. Thus, pivoting results in an 
(approximate) linear dependence between the logarithm of 
the flux (not flux itself) and the spectral index, as illustrated 
by the dashed lines in Fig. |7| 

Al Moments of the flux 

If pivoting variability occurs on time scales shorther than a 
given observation, the measured average spectrum will con- 
sist of the initial spectrum (in particular, a power law itself) 
times the average of /, /. If Ep = constant and the distri- 
bution in time of 5 is uniform from — Ar to -|-Ar, the flux 
average is given by, 



sinh X 



o 



(A2) 



where 



x = Arln{E/Ep). (A3) 

If the distribution of 5 is Gaussian, i.e., oc exp[— ((5/Ar)^/2] 
(where now Ar is the standard deviation of the distribution 
of S), we have, 

/ = e^'/' = l + y +0(x*). (A4) 

Note that the concave form of these spectra yield strong 
departures from the initial power law ai E <^ Ep. Pivoting 
on short time scales can account (at least in part) for the soft 
X-ray excesses commonly seen in Seyferts and BH binaries. 
Pivoting may also explain an apparently concave part of 
the extragalactic 7-ray spectrum (Stecker & Salamon 1996), 
although it requires the pivot energy to be within that part 
(which issue was not considered by those authors). 

This variability will also contribute to the flux variance, 
= {(,/ — .f)^). The variance normalized to the average flux 
is given by 

^ = a;cotha; - 1 = hO{x ), (A5) 

3 

for the uniform distribution of 5, and 
P 

for the Gaussian distribution of 5. With addition of some 
energy-independent variability, ctqi the ratio {a + o-q)// 
reproduces well the fractional variability of Cyg X-1 in the 
hard state on long time scales, see Fig.|5] 

Then, the skewness, s — {(/ — f)^) ja^ , is given by, 

6 -f x + 3 (x coth a:: — 3) coth x\ 



a 
72 



\ = x'' ^0{x^), 



(A6) 



3 (t coth a; - 1)^''^ 
^+0(.^ 



for the uniform distribution of 5, and 

1/2 



(A7) 



(A8) 
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for the Gaussian distribution of 5. This quantity may be use- 
ful for testing whether a given variability pattern is related 
to pivoting. Also, the kurtosis, S — {(/ — ^ 3, is 

given by, 



(e^%2) 



- 6 = 16x^ + 0(a 



(A9) 



for the Gaussian distribution of 5. 

On the other hand, it is rather unlikely that the pivot 
energy is completely constant in an astrophysical system. 
Instead, a range of pivoting energy is expected from some 
physical constraints on the emitting plasma, see fig. 3 in 
ZGOl and fig. 14 in Z02. Therefore, we also consider a 
case when the pivot energy is distributed log-normally, 
i.e., given by ln_Ep — ln_Ep — e, where e is distributed 
oc exp[— (e/Ap)^/2] and Ap is the standard deviation of this 
distribution. When S is also distributed normally, the aver- 
age departure from the original spectrum is given by. 



g:rV[2(l-A2)l 
(1- A2)l/2 
1 

- (1_ A2)l/2 

where 

A = ArAp < 1, 



2(1 - A2)3/^ 



x = Ar HE/E°). 



The variance is then. 



2 e' 
a — 



2X-' /{l~4,A^) 



(l_4A2)i/2 (l-2A2)i/2' 
and the variance normalized to the average flux is. 



(AlO) 



(All) 



(A12) 



A^ 



exp 



(1-f 2A^ 



(1 -4A2)(1 - A2) 



(1 - 2A 



2^1/2 



exp 



(1-2A2)(1- A2 



(1-4A 



2\-l/2 



1 -2A 



2\-l/2 



+ 



1 + 2A2 



(l-4A2)3/2 (i_2A2)3''2 



x^ + Oix"), (A13) 



for A < 1/2. Note that a finite range of the pivot energy 
results in a finite variance at i? = Ep, as given by the third 
line of equation IIA13t , in contrast to the case with the fixed 
pivot, equation (IA6II . Thus, the above dependence fits well 
the long-term behaviour of Cyg X-1 in the hard state with- 
out an additional constant, see Fig. |^ 
The skewness is given by. 



{(/ - If) = sa' 

9a;^/(2-18A2) 



+ 



2e' 



,3i^/(2-6A2) 



3e^ 



5a;^/{2-10A^) 



(l-9A2)l/2 (l-3A2)l/2 (l-5A2)l/2 



(A14) 



for A < 1/3. 

Above, we assumed that variations of the pivot energy 
are uncorrelated with those of the spectral index. Note that 
a correlation between these quantities would result in time 
lags between light curves measured at different energies (Ko- 
tov, Churazov & Gilfanov 2001). 



A2 Moments of the logarithm of flux 

The moments of the distributions above assume especially 
simple forms for the logarithm of flux. We define, 



g = lnf ^5lniE/Ep) 



(A15) 



where / is a departure from the initial spectrum, as above. 
Then, for a Gaussian distribution of 5 and a constant Ep, 
the distribution of g is Gaussian itself, with 



9 - 



0, 



0, 5 = 0. 



(A16) 



Thus, the apparently non-Gaussian distribution of the lin- 
ear fiux (equations |A8) - |X9] 'l becomes completely Gaussian 
when its logarithm is used. Such log-normal distributions 
are, in fact, very common in natural phenomena ranging 
from terrestrial lightning to 7-ray bursts (McBreen et al. 
1994). Note then that a positive skewness of the distribution 
of a linear variable (e.g., Leighly 1999) is not necessary an 
indication of non-Gaussianity of the underlying processes. 

In the case of Gaussian distributions of both S and e, 
g is a product of two variables {5 and e -I- constant), each 
having a Gaussian distribution, and 



5 = 0, 



0, s- 



6A2(A^ +2a;^ 
(A2 2^2)2 



.(A17) 



Note that a and S are not simply products of the corre- 
sponding quantities for each of the two variables due to a 
non-zero center value of the latter. Then, the resulting dis- 
tribution of g is not completely Gaussian, and, in particular, 
S>0. 
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